The PAD-US-AR dataset: Measuring accessible and recreational parks in the contiguous United States

Most spatial epidemiological studies of nature-health relationships use generalized greenspace measures. For instance, coarse-resolution spatial data containing normalized difference vegetative index (NDVI) values are prominent despite criticisms, such as the inability to restrain exposure estimates to public and private land. Non-threatening natural landscapes can improve health by building capacities for health-promoting behaviors. Recreational and accessible parks may best activate such behaviors. We curated the Parks and Protected Areas Database of the U.S. (PAD-US) to identify parks that are accessible for outdoor recreation. Our title adds “AR” to “PAD-US” where A = Accessible and R = Recreational. We validated the PAD-US-AR by comparisons with greenspace datasets and sociodemographics, which demonstrated its uniqueness from other commonly employed metrics of nature exposure. The PAD-US-AR presents reliable estimates of parks in the contiguous U.S. that are accessible for outdoor recreation. It has strong associations with home prices, shares of female residents, and shares of older residents. This dataset can accompany other nature exposure metrics in environmental epidemiology and allied research fields.


Methods
We curated the PAD-US-AR 48 dataset from the USGS Protected Areas Database of the U.S. V2.1 (PAD-US V2.1) 47 . The PAD-US is published by U.S. Geological Service in collaboration with Boise State University and through coordination with Federal, State, and non-governmental organizations that provide and verify the data. Its original release was in April 2009. Updates were made in 2010,2011,2012,2016,2018, 2020, and July 2022. Data on the completeness of the V2.0 dataset, which occurred before the V2.1 dataset used here, are available at https://www.protectedlands.net/ frequently-asked-questions-about-pad-us/. In brief, 14 states had over 95% coverage of parks and protected areas, 26 states had 80-95% coverage, and the remaining 8 states in the contiguous U.S. had <80% coverage. Updated coverage statistics for V2.1 are currently unavailable.
The PAD-US is a regularly updated geographic information system (GIS) spatial dataset that compiles the best available data provided by U.S.-based land management agencies and organizations. It strives to be a complete inventory of public land and other protected areas in the U.S. Included areas are those preserved for biological diversity and other natural, recreation, historical, or cultural uses and managed for these purposes through legal or other effective means 47 . Some areas consist of small land parcels with building footprints that occupy most of the area. These are not readily identified with the PAD-US V2.1 metadata. The location designation field (Loc_Ds) offers some clues with values such as "cultural arts center" and "National Register of Historic Places. " The number of unique values (N = 1,675) in the designation, easement, and fee areas of the PAD-US V2.1 limits precise identifications and removal of such areas.
The PAD-US V2.1 release became available in September 2020 and included notable updates from previous versions. These included integration of the TPL ParkServe dataset, Census American Indian/Alaskan Native Areas, Ducks Unlimited protected areas, and federal land ownership updates, among others. The PAD-US V3.0 was released in early July 2022 and contained minor updates that we expected to influence the curation process of the PAD-US-AR very little. For a complete description of version updates, see https://www.usgs.gov/ programs/gap-analysis-project/pad-us-data-history.
The PAD-US has been used for conservation mapping [49][50][51][52][53][54][55] and noise research 56,57 . These studies have identified that Western U.S. National Monuments provided jobs and economic growth after establishment 52 , counties with greater coverage of protected areas with strict conservation status (i.e., Wilderness Areas and National Parks) are associated with higher average noise levels 56 , and anthropogenic noise is common in many U.S. parks and protected areas 57 . We are also aware of a few nature-health studies that have utilized the complete PAD-US dataset 58,59 . In studies by Tsai and colleagues, the authors identified park locations and ground-truthed results with Google Maps and county/municipal data to identify park entrances. The PAD-US was used to calculate descriptive sample characteristics or covariates in models with other measures of nature exposure (i.e., tree cover and greenery), so associations between health and the PAD-US were not reported.
The opportunities and lack of precedent for curations of the PAD-US prompted us to define which types of parks and protected areas in the dataset were both accessible and recreation-oriented. Based on discussions We excluded the following designations (see the paragraphs below for rationales): 1. Department of Energy, Department of Defense, and Bureau of Reclamation lands 2. Marine areas managed as Marine Protected Areas by the National Oceanic and Atmospheric Administration, or Bureau of Ocean Energy Management, among others 3. Proclamation areas, which are boundaries of national lands used for administrative purposes that overlap with large areas of public lands that are not all available to the public 4. Fish hatcheries and other lands used for water rights with regulated hunting 5. National Park easements (i.e., lands paralleling but not including the Appalachian Trail and not used by the public) 6. Joint management areas (i.e., university research stations) 7. Non-governmental organization lands (aside from conservation easements) 8. State trust/land survey lands 9. American Indian Lands 10. Other areas with unknown access or closed public access (i.e., limited to coordinated programs and research) Restricting the PAD-US to these categories was a sequential process starting with the four terrestrial PAD-US domains (Fig. 1  www.nature.com/scientificdata www.nature.com/scientificdata/ and proclamations (boundaries of administrative areas). For further information on these domains, see http:// www.protectedlands.net/pad-us-technical-how-tos/.
Our first step was to exclude all proclamation lands in the PAD-US. These administrative boundaries are not ownership lines but are instead the outer boundaries of areas used by land managers for planning regardless of internal ownership. They could but will not necessarily be publicly managed in the future. Some commercial mapping providers (i.e., Google Maps, Esri USA Parks) incorrectly use these boundaries to show protected areas and, in doing so, often show large areas of private lands as part of public lands.
Next, we excluded lands described as closed to public access in the PAD-US. Alternative classifications include open to public access, restricted, which denotes a permit may be needed, or unknown. We temporarily retained unknown access areas for further consideration since large areas of the intermountain west are designated as such. For example, the Great Salt Lake, UT, is the state's largest water body and a recreation destination for boating, swimming, and sunbathing.
The subsequent step was refining lands labeled as unknown access in the PAD-US. Decisions were made based on the assigned land manager. City lands (Code = CITY) were included since many greenways were under this classification. County lands (CNTY), which described nearly 250 polygons run by the City of New York for parks and recreation in the city and upstate, were included. Similarly, regional agency land (REG) covered over 400 polygons concentrated in Chicago and Los Angeles suburbs used for parks and recreation; these lands were retained. State Department of Conservation (SDC) and State Department of Natural Resource (SDNR) lands were also retained. These included over 5,000 polygons across the country, including the Great Swamp Management Area, RI, an important area for birding and open to the public, and the Great Salt Lake. State Department of Land (SDOL) areas were retained, as they included approximately 30 polygons used by the public for hiking in Northwestern states. State Fish and Wildlife (SFW) lands included urban areas with trails along waterways and were retained. State Parks and Recreation (SPR) lands were retained and covered public recreational areas in Maine. Tennessee Valley Authority (TVA) and Army Corps of Engineers (USACE) areas covered large reservoirs with important water-based recreation resources and were retained. The presence of such waterbodies, which provide public recreation to millions of visitors annually 61 , required us to retain the entire census geographies despite evidence that removing areas covered by water can lead to more precise and realistic sociodemographic analyses 62 . Last, U.S. Forest Service (USFS) lands were retained as they included several recreational areas in Virginia.
All other areas with unknown public access in the PAD-US were deemed not accessible to the public or used for public recreation and therefore excluded. This conservative approach reduced the chances of misclassifying large tracts of land that were likely inaccessible. For example, Department of Defense (DOD) lands included ammunition plants, Department of Energy (DOE) lands included nuclear test sites, and National Oceanic and Atmospheric Administration (NOAA) lands included estuarine research reserves. Non-governmental organization (NGO) lands included nearly 17,500 polygons in the Rocky Mountains but covered too many conservation types to determine whether these were open to the public. American Indian Lands (TRIB) were on reservations and could not be assumed to be accessible and used by the public.
The final step in curating the PAD-US-AR dataset was determining how to approach the polygons in the Western and Midwestern states that were left over from the Public Land Survey System (designation = SRMA). Most of these lands follow a grid pattern and are not used for outdoor recreation. However, some state trust lands include important parks, such as DuPont State Forest, NC, a popular destination for mountain biking, hiking, swimming, and visiting waterfalls. Three of the authors conducted online searches of the uses of these lands using online resources (i.e., State Department of Natural Resource portals) for each state and selected which to include or exclude. The corresponding author also discussed these decisions with three outdoor recreation professionals living in the western U.S. Based on this examination, we removed state trust lands from Arizona, Colorado, Idaho, Louisiana, Mississippi, Montana, Nevada, New Mexico, North Dakota, Oklahoma, Oregon, South Dakota, Texas, Utah, Washington, and Wyoming.
To obtain census tract and county exposure estimates, we calculated the percentage of the PAD-US-AR covering each geographic unit. Tract-level estimates included a 0.5-mile buffer around each tract to acknowledge the opportunities for park access for residents living near the tract boundaries. Similar thresholds have been used in past research [63][64][65] and are recommended as U.S. park access standards by several nonprofits (e.g., Trust for Public Land, www.10minutewalk.org). This threshold is primarily used in urban areas and may be most relevant to those areas where most people live and where tract sizes are smaller.
No buffering was conducted around counties. Counties are >300% larger than tracts, on average. In our dataset, the median county area was 1,614 km 2 , while the median tract area was 5 km 2 . Counties are also jurisdictions of local governments, whereas tracts do not represent any administrative boundaries. For these reasons, we avoided buffering counties, which often have parks and recreation departments managing parks within their borders. www.nature.com/scientificdata www.nature.com/scientificdata/ Geopackage and shapefiles include vector polygons with the original metadata from the PAD-US V2.1. For a complete listing of variables, please visit https://www.usgs.gov/programs/gap-analysis-project/ pad-us-data-manual. In brief, the data include the name of the parcel; feature class (in the PAD-US-AR, the options are designation, easement, or fee); type and name of management agency (i.e., federal, state, American Indian Lands, or local government); designation (i.e., conversation easement vs. National Park); conservation protection level as designated by the International Union for the Conservation of Nature (IUCN); state name; and geographic size.

Data Records
Spreadsheets include geographic identifiers (i.e., FIPS codes or GEOID) and percent park cover. These are provided as Microsoft Excel files (.xlsx) and text files (.txt) to maintain leading zeros in the geographic identifiers. Park cover ranges from 0 (no parks) to 100 (complete park cover). Tract estimates are provided for park cover within the boundaries of each tract and the 0.5-mile buffered tract boundaries.

technical Validation
The PAD-US-AR 48 dataset presents park cover from nearly 250,000 spatial units and 1,900,000 km 2 in area across the contiguous U.S ( Table 2). Histograms of the data within counties and tracts and by census region are presented in Figure S2. Distributions were right skewed in all regions except Northeastern and Western counties. Northeastern counties showed a flat distribution until approximately 20% cover. Higher levels of cover were present in few counties. Western counties showed a roughly flat distribution until around 80% cover, after which the number of counties with higher cover levels was small.
Comparisons with the source dataset (PAD-US V2.1) are available for each census region in Figs. 2-5. Large areas of Maine, southeast Pennsylvania, central/western Massachusetts, and northern New Hampshire were excluded from the PAD-US-AR because they were private conservation easements, watersheds with closed access as listed in the PAD-US V2.1, or otherwise unknown public access. Swaths of the Dakotas were removed as conservation easements used for wildlife management with uncertain public access. Lands in Oklahoma arranged on a gridwork were removed as state school lands typically leased out for agriculture and mineral resource purposes. A gridwork of land parcels in Montana, Wyoming, Colorado, Arizona, and New Mexico was removed as state trust lands managed for timber, surface, and mineral resource extraction. Similarly, larger parcels of state trust lands in Western Texas were excluded. Other large parcels of lands excluded were over 560,000 acres in central Idaho, 860,00 acres in southern Nevada, and nearly 200,000 acres in southern South Carolina managed by the Department of Energy; approximately 550,000 acres at Vermejo Park Ranch managed by Ted Turner Reserves, Inc., and 133,000 areas of the Stronghold District of Badlands National Park in western South Dakota owned by the Oglala Sioux Tribe under agreement by the National Park Service.
Next, we compare the PAD-US-AR dataset with other park datasets, nature exposure metrics, and sociodemographic characteristics. The value of comparing the PAD-US-AR with other park datasets is to determine whether the PAD-US-AR differs from already available datasets. Park dataset comparisons were made by tallying the number of geographic polygon units and calculating the total cover after dissolving all polygon units (to account for some polygons overlapping each other) in census regions.
The value of comparing the PAD-US-AR to nature exposure metrics is to determine whether park cover differs from other standard exposure estimates. We employed two measures of NDVI (annual averages and summertime highs) and tree canopy cover, which were derived from raster images and averaged across geographic units (tracts or counties). NDVI values were retrieved and processed in Google Earth Engine (GEE) using cumulative annuals or summertime highs (June-August) from 250 × 250 m 16-day MODIS images averaged over five years (2015-2020) after extracting cloud cover and water pixels. Tree canopy data were retrieved from the 2019 National Land Cover Database (NLCD) release 66 , which provides cover estimates ranging from 0 to 100% for each 30 × 30 m pixel in 2016. This release was the most recent available during data retrieval (September 2022). To identify whether the PAD-US-AR was unique from these other estimates of nature exposure, we examined bivariate correlations between each metric and the PAD-US-AR.  www.nature.com/scientificdata www.nature.com/scientificdata/ Last, we examined sociodemographic correlates of park cover measured through the PAD-US-AR to inform what confounding factors should be considered when modeling associations between park cover and human health. Sociodemographic characteristics were retrieved from 2015-2019 American Community Survey (ACS) estimates from the U.S. Census at the county and tract level 67 . We selected 14 variables (Table S2) based on existing literature examining correlates of greenspace, especially in studies focusing on socioeconomic and racial disparities in access to these spaces [68][69][70][71][72] . Attempts at incorporating median household income alongside other measures resulted in multicollinearity, so this variable was excluded from the primary analyses but considered in a sensitivity analysis. We examined the results of generalized linear mixed models (GLMMs) with gamma distributions and U.S. states as random effects to account for the non-normal distribution of the outcome variable and the hierarchical nature of the data (counties and tracts within states). Models were run with complete data for 100% of counties (N = 3,108) and 97.3% of tracts (N = 70,580) in the contiguous U.S. circa 2019.
Stratified analyses using more urbanized counties (≥50 people/km 2 ) and tracts (≥1,000 people/km 2 ) were conducted to compare results with past research and inform future scholarship with the PAD-US-AR. There is no consensus on differentiating more vs. less urban areas in nature-health research 19 . Between 1,000 and 1,999 people/km 2 is a common cut point 19 . We attempted to apply that cut point to both units of analysis (tracts and counties), which split the number of tracts roughly in half (n = 32,929 as more urban). In contrast, this cut point resulted in too few counties to conduct sufficiently powered analyses (N = 45 as more urban). We attempted the 300 people/km 2 cut point recommended by the European Union 73 and used in a recent U.S. study on the association between park cover, park use, and mental health 74 . This continued to produce small sample sizes: N = 43 for the Northeast, 30 for the Midwest, 93 for the South, and 16 for the West. A cut point of 50 people/km 2 produced reasonable sample sizes for most regions (N = 121 for the Northeast, 178 for the Midwest, 386 for the South, and 58 for the West). Applying this 50 people/km 2 cut point to counties also produced maps that approximated the location of the Census classification of urbanized areas (https://www.census.gov/programs-surveys/geography/ guidance/geo-areas/urban-rural.html; Figure S2). This urbanized area classification scheme has been used to create other datasets on environmental exposure estimates, such as urban heat island vulnerability 75 . GLMMs were used in these stratified analyses except in the Midwest, where standard linear regression models were run to avoid singularity resulting from few urban counties per state in the random effect term.

Comparison of the PAD-US-AR percent park cover dataset to other park datasets.
Descriptive statistics for each park dataset are provided in Table 2, and maps of park cover are provided in Figure S3. The PAD-US-AR 48 covers 51.6% of the acreage in the PAD-US V2.1 dataset. The PAD-US-AR acreage is larger than the acreage of USA Parks and ParkServe but smaller than the OSM datasets when leisure and boundary tags  www.nature.com/scientificdata www.nature.com/scientificdata/ are combined. Bureau of Land Management (BLM) lands are mainly absent from the USA Parks and ParkServe datasets but are partially included in the OSM datasets and prominent in the PAD-US-AR. This is particularly noticeable in Nevada, western Utah, and Wyoming. These areas include such popular recreation attractions as the Grand Staircase-Escalante National Monument, UT, and the Grand Canyon Parashant National Monument, AZ. These collectively encompass nearly 3,000,000 acres (around twice the size of Delaware), attract more than 150,000 visitors annually for hiking, backpacking, and camping, and have received thousands of 5-star reviews on Google Maps. This high number of reviews shows their popularity and visibility in the public sphere. Other notable areas include off-highway vehicle (OHV) trails, such as the Little Sahara OHV Area, UT, which offers driving/riding on a 700-foot drivable sand dune, 30,000 annual visitors, four campgrounds, and approximately 62,000 acres. Most popular mountain biking and OHV riding trails around Moab, UT (except for the Slick Rock Trail System) are also BLM lands excluded or with limited coverage in datasets other than the PAD-US and PAD-US-AR. These results demonstrate that the PAD-US-AR presents a selected sample of the PAD-US dataset with differing coverage from pre-existing park cover datasets.

Comparison of the PAD-US-AR to other nature exposure measures. Descriptive statistics for park
cover and other nature exposure metrics are presented in Table 3. Maps of each metric are provided in Figure S4. Distributions of nature exposure metrics are available in Figures S5-S7.
Associations between the PAD-US-AR 48 and NDVI varied across geographies and seasons (Fig. 6). Park cover was negatively associated with NDVI at the county level. Pearson correlation coefficients (r[95% confi-     www.nature.com/scientificdata www.nature.com/scientificdata/ These findings demonstrate that the PAD-US-AR 48 presents a unique exposure estimate from metrics of nature exposure. Plant-rich landscapes, or "greenspaces, " do not capture all aspects of open recreational spaces and nature-rich landscapes 10 . Correlations between nature exposure metrics vary in size and direction based on the unit of analysis (counties vs. tracts) and geography (regions of the country and nationwide analyses).

Comparison of the PAD-US-AR to sociodemographic characteristics.
A listing of the sociodemographic characteristics considered in analyses is provided in Table S1. Descriptive statistics for each variable are presented in Tables S2-S6. Maps of the distribution of these variables are provided in Figure S8. Multivariate associations between the PAD-US-AR 48 and sociodemographic characteristics are shown in Fig. 7 and Table S7. These results were derived from GLMMs with gamma distributions and U.S. states as random effects to account for the non-normal distribution of the outcome variable with minimal multicollinearity (Table S8).
Park cover was more strongly associated with sociodemographic characteristics at the county level than at the tract level. Around 30% of the variance in countywide park cover was explained in U.S. regions after accounting for state random effects (conditional R 2 Northeast = 0.29, R 2 Midwest = 0.31, R 2 South = 0.23, R 2 West = 0.38). Variance explained within counties across the country was over 60% (R 2 Nationwide = 0.63). Variance explained at the tract level was closer to 10%-20% (R 2 Nationwide = 0.19, R 2 Northeast = 0.09, R 2 Midwest = 0.08, R 2 South = 0.12, R 2 West = 0.18).

Fig. 7
Regressing sociodemographic characteristics on the PAD-US-AR park cover dataset within counties and tracts. Notes: GLMMs with gamma distributions and U.S. state random effects, except for Midwestern counties, which report standard linear regression results. Standardized betas and 95% confidence intervals are shown. Differing symbols represent statistical significance (p-value): an empty circle is shown for p > 0.05, a filled-in circle for p < 0.05; a triangle for p < 0.01; a square for p < 0.001. X-axes are on different scales. Sensitivity models with median household income substituted for other socioeconomic variables are provided in Figure S9. www.nature.com/scientificdata www.nature.com/scientificdata/ Three sociodemographic characteristics showed fairly consistent associations with park cover. On average, areas with greater shares of older adults (≥65 yrs) had more park cover. Areas with higher median home values also had more park cover, except in the Northeast. Last, areas with greater shares of female residents had less park cover on average, except in Northeastern and Southern counties. Two other sociodemographic characteristics showed consistent associations within either county or tract samples. First, counties with lower Gini index values (lower inequality) had more park cover on average. Secondly, tracts with higher unemployment rates had more park cover on average.
Associations between the PAD-US-AR and other sociodemographic characteristics varied by region. Park cover in Northeastern counties was concentrated in areas with lower rates of income inequality, high school graduation, and natural resource employment. Park cover in Midwestern counties was greater in areas with higher poverty rates, income inequality, and unemployment. Park cover in Southern counties was higher in areas with greater population densities or higher rates of poverty and lower rates of income inequality, natural resource employment, or non-Hispanic Black residents. Western counties showed greater park cover in areas with more poverty, higher shares of college degree holders, less income inequality, and lower shares of non-Hispanic Asian residents. Within tracts, park cover was higher in densely populated Northeastern areas but lower in densely populated areas throughout the rest of the country. Tract-level park cover was higher in areas with greater shares of residents employed in natural resource professions in the West and Northeast, while the opposite was found in the South and Midwest; in these areas, park cover was lower in areas where greater shares of people worked in natural resources professions. Park cover was higher in Midwestern and Southern tracts with greater shares of non-Hispanic Asians, whereas park cover was lower in Western and Northeastern tracts with greater shares of non-Hispanic Asians. In summary, park cover was associated with many sociodemographic characteristics, but the strength and direction varied by geography and unit of analysis.
Multivariate associations between the PAD-US-AR and sociodemographic characteristics in urban areas are presented in Table S10. In most cases, median home value continued to show strong positive associations with park cover. One exception was observed in Midwestern tracts, where median home value was negatively associated with park cover. Percent female no longer predicted park cover except in Southern tracts. Shares of older adults also predicted park cover in only a few urban cases; significant positive associations were observed only in nationwide and Northeastern tracts. Percent Non-Hispanic Asian residents emerged as a predictor in several models, but the direction of the associations differed. Nationwide models showed negative associations, while Midwestern counties and tracts and Southern tracts showed positive associations. County-level models of urban areas continued to predict the variance explained by park cover better than tract-level models of urban areas. Alternative models substituting median household income for other socio-economic indicators found mixed relationships between this variable and park cover (Table S9, Figure S9).

Usage Notes
We present a new potential indicator of outdoor nature exposure for the contiguous U.S: the location of parks intended to be accessible for recreation. This dataset allows researchers to examine the number of outdoor recreation areas meant for public use around geographic units of interest (i.e., homes, neighborhoods, and transit routes). Other commonly-used metrics -like moderate/coarse resolution NDVI and tree canopy cover datasets -cannot identify whether the areas are managed for public recreational use. The PAD-US-AR 48 is unique from these other metrics, as determined by the correlations presented above (Fig. 6).
The PAD-US-AR also differs in coverage from pre-existing park cover datasets. These differences were observed when tallying the geographic polygon units and calculating the total cover after dissolving all polygon units to account for some overlapping units. The reasons to utilize the PAD-US-AR dataset rather than these other options include the PAD-US-AR source data (PAD-US V2.1) being validated by the agencies managing the land, our systematic examination of what is accessible for recreation, and the clarity and transparency in its curation. The potential for park cover to not match park access for all residents in a county or tract remains high, as in any area-level exposure estimate [80][81][82] . Individual-level estimates should be calculated from the boundaries or centroids of park polygons along road or pedestrian networks when geolocated data for homes, schools, workplaces, or activity spaces are available.
The chances for residual confounding in area-level studies with the PAD-US-AR dataset exist if multivariate models do not control for sociodemographic characteristics of the areas encompassing parks. The PAD-US-AR has the most robust associations with home prices, shares of female residents, and shares of older presents. These should be statistically controlled in models using the PAD-US-AR as an independent variable or covariate. Other measures of socioeconomic status (i.e., median household income) might be insufficient to avoid residual confounding in ecological studies with PAD-US-AR data.
Since the PAD-US-AR was curated nationwide, it is most appropriate for use at larger geographic scales (i.e., regional and national). Studies focusing on smaller geographic contexts, such as within individual cities or states, should partner with local land management agencies and recreation departments to ensure PAD-US-AR data accurately represent all parks and protected areas managed for public outdoor recreation. Since ownership boundaries and land acquisitions can change annually, local land management agencies might also be able to identify new parks that aren't present in the PAD-US-AR. Smaller-scale analyses may allow manual selection of land parcels with building footprints that occupy most of the area.
The PAD-US-AR may be best conceived as the minimum park coverage level. We excluded the approximately 35,000 areas covering over 42,000 km 2 with unknown public access in the PAD-US. Some private parks, such as golf courses or community parks restricted to residents who pay homeowner association fees, can provide opportunities for outdoor recreation that activate the same health-promoting pathways as public parks.
www.nature.com/scientificdata www.nature.com/scientificdata/ People living in the counties and tracts presented in the datasets may have more access to outdoor recreational opportunities than suggested by the PAD-US-AR.
As the nature-health literature expands, exposure estimates are expected to develop and be refined. The PAD-US-AR presents a significant advancement in this body of literature by offering researchers an assessment of where parks are available for outdoor recreation.

Code availability
R (4.1.2) was used to generate the PAD-US-AR 48 dataset and results. QGIS (3.18.3) was used to create the maps. Scripts and source data to reproduce results are available on OSF (https://osf.io/pwdsg/).